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Prefiice 

This  thesis  on  lov/  thrust  orbital  transfer  v>/as  su<j<jested  by  an 
article  which  appeared  in  a IS7I  issue  of  The  Journal  of  the  Astronau- 
tical  Sciences  by  J.L.  Starr  and  R.O.  Sugar.  While  their  analysis  was 
limited  to  transfer  between  coplanar  orbits,  It  was  rny  desire  to  extend 
this  problem  to  the  more  general  case  of  transfer  between  orbits  of  dif- 
fering inclination. 

The  application  of  modern  control  theory  to  such  a non-linear  system 
ultimately  necessitated  cor’puter  utilization.  As  is  often  the  case, 
this  computational  phase  of  engineering  was  not  only  t Ime-consurning,  but 
frustrating  as  well. 

In  my  case,  the  latter  difficulty  was  eventually  surmounted  through 
the  encouragefflent  and  understanding  of  my  v;lfe,  Glenda.  The  knowledge 
and  advice  of  my  advisor,  Professor  Gerald  M.  Anderson  of  the  Air  Force 
Institute  of  Technology,  was  also  Invaluable  towards  the  successful  com- 
pletion of  this  work. 
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Abstract 

Minimum  fuel  continuous  low  thrust  orbit  transfer  is  an  optimal 
control  problem  v/hlch  requires  the  application  of  iterative  numerical 
methods  to  solve  the  resulting  two  point  boundary  value  problem.  A 
combination  gradient-neighboring  extrctnal  alnorithm  gives  accurate  re- 
sults when  both  the  initial  and  terntnal  orbits  are  coplanar.  These 
results  are  then  used  to  extend  the  problem  to  include  out-of-plane 
transfers  between  orbits  of  V'.'idely  varying  inclination.  Control  his- 
tories, in  terms  of  thrust  direction,  are  plotted  as  a function  of 
time  for  both  coplanar  and  non-coolanar  transfers.  Steering  programs 
for  coplanar  transfers  are  characterized  by  a rapid  thrust  reversal 
which  is  dependent  on  eccentricity  for  the  time  of  its  occurrence. 

In  the  latter  case,  out-of-plane  tlirusting  is  primorily  a function  of 
the  inclination  change  between  orbits.  On  the  basis  of  these  results, 
the  particular  numerical  approach  used  provides  an  effective  method 
for  generating  optimal  low  thrust  orbital  transfer  guidance, 
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MINIMUM  FUEL  CONTINUOUS  LOW  THRUST  ORBITAL  TRANSFER 
* • Introdiict  Ion 

While  our  present  chemical  prooulslon  systems  have  led  to  tremen 
dous  accomp I I shnents  In  the  field  of  space  exploration,  it  seems  ap” 
parent  that  the  day  of  large,  bulky,  impulsive-thrust  rockets  will 
soon  be  ending.  The  ever  increasing  costs  of  these  systems  coupled 
with  the  depletion  of  our  energy  resources  will  soon  make  the  develop 
ment  of  a low  thrust  propulsion  system,  such  as  an  ion  or  nuclear  en- 
gine, imperative. 

The  purpose  of  this  thesis  is  to  examine  the  minimum  fuel  orbit 
transfer  problem  In  the  context  of  a continuous  low  thrust  propulsion 
system.  Since  the  propulsion  system  is  assumed  to  give  constant 
magnitude  low  thrust,  a minimum  fuel  transfer  becomes  synonymous  with 
a minimum  time  transfer.  While  this  particular  problem  has  been 
analyzed  extensively  in  terms  of  coolanar  orbital  transfer  (Ref 
204),  the  overall  purpose  of  obtaining  numerical  solutions  here  Is 
tv/o-fold:  investigate  the  possibility  of  using  a combined  gradient 

a''  ' neighboring  extremal  algorith:.'.  for  coplanar  transfer,  and  apply 
th  -.oplanar  results  as  a reference  input  in  extending  the.  problem 
to  non-coplanar  transfer. 

The  problem  is  to  determine  the  optimal  thrust  direction  to 
transfer  a vehicle  from  some  point  on  an  initial  orbit,  to  some  de- 
sired terminal  orbit  in  minimum  time,  using  a constant  magnitude  lov/ 
thrust  propulsion  system. 

In  Chapter  II  both  the  initial  and  terminal  orbits  are  assumed 
to  be  coplanar  .ind  coapsidal,  while  the  problem  is  extended  to  a non- 
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coplanar  terminal  orbit  In  Chapter  III.  In  each  case,  the  equations 
of  motion  are  developed  from  thf  restricted  two-body  problem  with  addi- 
tional terms  to  account  for  vehicle  thrust.  For  this  analysis  a 5000 
pound  vehicle  with  a thrust-to-welght  ratio  of  0.0125  Is  assumed. 

In  this  paper  the  problem  is  a.5oroached  with  optimal  control 
theory.  The  thrust  direction  becomes  the  control  variable  and,  through 
the  equations  of  motion,  a Two  Point  Boundary  Value  Problem  (TP15VP)  Is 
formulated.  Numerical  solutions  are  then  obtained  to  the  resulting 
TP3VP  using  a combination  gradient-neighboring  extremal  algorithm  for 
the  case  of  coplanar  orbits.  Solutions  for  the  non-coplanar  orbits 
are  obtained  by  a neighboring  extremal  algorithm  using  the  coplanar 
results  as  an  initial  reference. 

The  results  are  presented  graphically  In  Chapter  IV  In  terms  of 
the  optimal  transfer  trajectory  and  optimal  control  history.  The  con- 
clusions are  discussed  in  Chapter  V. 
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Specifying  the  seml-latus  rectum  and  eccentricity  of  the  Initial 
and  final  orbits  as  p^,  e^,  pp  and  e^,  respectively,  and  the  true 
anomaly  of  departure  from  the  Initial  orbit  as  O^,  the  Initial  condi- 
tions are  given  by 


1 + cos  0^ 

Vr  -fo!klL!l®o 

,^0 


V0O  - 

•■o 


where  h^  Is  the  specific  angular  momentum  of  the  initial  orbit. 
Tfie  terminal  constraints  arc  then  determined  as 


f " ; 


I + Cj:  cos  Op 
Ophp  sin  Op 


VOf  - ~ 

^ •‘f 


where  hp  is  tlie  specific  angular  momentum  of  the  terminal  orbit  and 
Op,  Lhe  true  anomaly  at  arrival  to  the  terminal  orbit,  Is  unconstrained, 
Using  the  Mayer  formulation  of  the  cost  function  In  terms  of  a 
minimum  time  trajectory 


J « t, 


whlch  yields  the  following  for  the  Hamiltonian 


s 

(i* 
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„ . i,  Vr  ♦ Xj  * X,  “-f".  - a--  + I-t""-  s 

+ Xk  ,•  rliys^  + I_52L.u.-n 
- mt  y 

(12) 

Tlie 

costotc  equations  thus  become 

. VJ.  , V02  . Vr,^ 

^ ',r2  ^ ’ f2  ^ 

(13) 

« 

X2  " 0 

(HO 

h - - h +7-H 

(15) 

^ \ 2V0  , ^ Vr  , 

-'4  y *2  ■ -|T*  A3  + X4 

(16) 

v/I  th  the 

following  transversal i ty  conditions 

Xi  (tf)  “ vj  + V3  ■— i- 

(17) 

x,(t,) . - 

(1  •♦■  Cf  cos  Of)^ 

(18) 

X3(tf)  * V2 

(19) 

A4(tf)  •»  V3 

(20) 

H(tf)  - - 1 

(21) 

whore  vi 

, V2.  and  V3  are  constant  Layrango  nul tipi lers. 

Mlnlmliiing  the  Hamlltoniiin  vvJtli  respect  to  u ami  enforcing  the 
strengthened  Lcgendre-Clchsch  condition  gives 
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sin  u* 


and  cos  u”  * 


/X32  + X|^2 


as  the  optimal  control. 


Summary  of  the  Coplanar  TPBVP 

Defining  the  state  vector  components  r,  0,  Vr,  and  VO  as  Xj , x^, 
X3,  and  Xj,,  respectively,  the  coplanar  TPEVP  becomes 


xi(0)  « r^ 


X2(0)  - Oj 


“1  “i* 



(m^  - rnt)  + X^^ 


Vr^ 


-X3X4 


x4(o)  » 


(m_  “ nit)  A, 2 + >,2 


?ii_  X + 


X.  (tp)  “ V.  + V — (27) 
>^\r 


X^  Xj 


- V 
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Xj  “ 0 


’ X4 

^3  “ ‘^1  j<Y  ^4 


2xk 


\ - -_r-  X 

X^  2 Xj  3 


^3 

+ X. 

Jt  ^ 


-PfOp  sin  X - 
XjCtp)  - f r 2r 


0 + 


Be  COS 


Xjfl" 


“Cphp  COS  x^p 

Pf 


X3(tp)  - 


X.  (tr)  “ V 
4 ' 3 


l//|  - X 


Pf 


1 '‘if 


1 + Of  cos  X^f 


'^2  " ^3f 


ophp  sin  Xyp 


- - 0 


Pf 


i|/  ••  X “ -J—  • 0 


Q » ll(tp)  + 1 » 0 


8 constants  of  integration 
+3  unknown  multipliers  v. 


*>  12  unknowns 


In  It  III  1 conditions 
+3  terninfil  constrni  'ts 
1.^1  ^transynrsn  1 i ty_  cpiidjt  jon£ 
« 12  cciuatlons 


(28) 


(29) 

(30) 


(31) 


(32) 

(33) 


(3M 


Thus  In  theory  It  is  possible  to  solve  this  TPDVP.  However,  clue 
to  the  non-linearity  of  the  equations  ami  their  o/pllctt  cJopeniJcnce  on 
tinio  an  analytic  solution  cc}ulcl  not  be  found. 


Grad  lent  Algor  I thn 

The  gradient  method  proved  to  bo  an  cxccilicnt  means  of  obtaining 

Reproducfd  from 

bett  availabU  copy.  WW 


GA/MC/7'»-5 


approximate  solutions  to  this  TPDVP.  The  particular  approach  used 
here  Is  the  method  of  multipliers  suggested  by  Hesteucs  (llcf  3*'10). 

It  consists  of  reformulating  the  optimal  control  problem  hy  forming 
an  augmented  penalty  function  of  the  form 

Jg  « tf  + (Sjj  4-  $22  '^2^  * S33  '^3^)  ^ 

(35) 


where  both  quadratic  and  linear  terms  In  >i.‘  have  boon  added  to  £q  (11) 
to  enforce  the  termlf^al  constraints  (itef  ^:^l03).  The  S terms  In  Eq 
(35)  ane  weighting  terms  which  may  be  varied  as  necessary  according  to 
the  terminal  error.  The  multipliers  bj  arc  updated  from  the  rela- 
tionship 

bj  - S|  (36) 


In  such  a mat.ner  as  to  Increase  In  magnitude  as  the  terminal  condi- 
tions become  closer  to  being  satisfied. 

This  change  In  formulation  ollmIn.ites  the  terminal  constraints, 
Eqs  (d) , (S) , and  (lO),  from  the  problem,  The  only  odJItlonal  change 
Is  then  through  the  transverse  1 i ty  conditions,  Eqs  (17),  (I"),  (19), 
and  (20),  which  now  have  the  form 


II 


If 


I + 0|:  cos  X^p  ^ 


+ S33 


hf  hp 

.J_*N  _J. — + bi  + h3 
xip  y Xip2 


hp 


^If 


(37) 


a 
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Xj  (tf)  - Sji  ^Xjj: 


I + Cf  cos  X2f 


C ZA-l-iy N 

^ (l  + cos  ^ 


/ sin  X*,  V 

(X3f  - -f_f_ fl.) 

Pf 


X “Ofhf  cos  . ■ PfCf  sin  x^f 

f’f  0 + Of  cos  Xjf)^ 


(■  cos  Xjf  * 

) 

(33) 

i x,Uf) 

■ s (" 
22  ^ 

sin  X,rv 

LfA  + b 

Pf  y 2 

(30) 

[;  Xjj  ( t j: ) 

" ^33 

hf 

^If^ 

) + ^53 

(40) 

Tlio  qrodlcnt  algoritlin  Is  described  throupb  the  follov-zlnq  steps: 

(I)  With  some  initial  ostimeto  for  u is  n function  of 

time,  o.g.,  u,-gf(t)  ■ 0,  the  original  state  equations, 

(l)  thru  (4),  arc  Intcqrntcd  forward  In  tine  until 
dJ,  d^J- 

...Jl  * 0 an(j a.  > 0 v/hicli  determines  an  Initial 


final  time,  tf. 

(2)  Eqs  (37)  thru  (*t0)  arc  used  as  initial  conditions 
with  b|  « 0 on  the  first  Iteration  to  Intoqratc  the 
adjoint  equations,  (13)  thru  (16),  backward  in  time 
using  the  non-optimal  reference  trajectory  from 
step  1 . 

(3)  Simultaneously  with  step  2,  the  gradient,  lly, 

Is  evaluated  at  each  point  along  this  rofcronce 


GA/MC/7^*-5 


trajectory. 

(^)  Since  to  first  order,  <SJg  “ lly  6u,  a one  dimensional 
search  Is  made  on  the  parameter  a where 

iSu  ■ “ally,  a > 0 (1*1) 

to  decrease  the  aunmentud  cost  function, 

(5)  The  control  is  iiodatcd  through 

«now  “ "ref  + (^2) 

wiicrc  becomes  the  new  estimate  In  step  1. 

(6)  The  b|  are  determined  from  Eq  (3C)  and  added  to  those 
used  In  step  2.  This  Is  an  Important  part  of  this 
procedure  because  as  the  terminal  constraints  become 
closer  to  being  satisfied,  the  linear  terms  In  Eq 
(33)  must  become  dominant  to  ensure  continued  con- 
vergence. 

(7)  Stops  I thru  C are  repeated  .ns  necessary  until  the 
terminal  constraints  are  satisfied  to  yield  the 
optimal  control,  u’‘(t},  and  the  minimum  time  tr.?ns- 
fer  trajoctoiy. 

This  niothod  is  quite  straightforward  In  approach,  but  as  with 
most  gradient  schemes,  It  re<iulrcs  many  Itcrotlons  and  Is  slow  to 
converge  nc.ir  the  optimum.  However,  a technique  used  In  the  actual 
computer  program  was  very  effective  In  ovorcof'lng  this  problem  and 
will  be  discussed  in  Chapter  IV.  Moreover,  the  results  obtained 
served  as  an  c;<cc1Icnt  reference  Input  to  the  neighboring  extremal 
algor  I thm. 
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Neighboring  Extrcnul  ,Al_gorJ  thm 

A major  obstacle  In  using  a neighboring  extremal  algorithm  Is  the 
necessity  of  a good  initial  estimate.  This  difficulty  Is  overcome  by 
using  the  results  of  the  gradient  method. 

The  particular  neighboring  extremal  algorithm  used  here  consir.ts 
of  estimating  the  unspecified  terninal  conditions  of  the  TPQVP  and 
generating  a transition  matrix,  which  is  used  to  adjust  the  terninal 
conditions  so  as  to  satisfy  the  specified  initial  conditions  (Ref  I: 
213).  This  Iterative  procedure  Is  as  follows: 

(1)  Estimate  tlie  following  eight  parameters;  x^j.,  x^p 

^3f’  ^4f’  ’■f*  Initially  this  Is 

dona  using  the  results  from  tho  gradient  method. 

(2)  The  eight  parometot's  in  step  1 determine  Xjf* 

^3f’  \f’  '^1’  '^Z'  '^'3’ 

(3M. 

(3)  Tlie  state  and  costate  equations  are  integrated 
backv;ard  from  tf  to  give  x^q,  x^q,  X3Q,  and  Xi^q 

which  arc,  in  general,  different  from  tlie  speci- 
fied Initial  conditions. 

ik)  Through  small  perturbations  in  tho  eight  para- 
meters of  step  1,  a transition  matrix  of  the  form 

5 (x^Q.  X3q.  C't,  'P3.  n) 

0 (x^j;,  ^2f*  ^3f*  ^4f*  '^l*  ^2’  *^3’  ^f^ 


is  numerically  generated  such  that 
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6 X (0) 

dijr 

3 (^io»  ^20’  ^30’  \o’  '^1' 

il'jt  >I'3»  Si) 

5 X (tp) 
dv 

B 

3 (^lf»  ^3f*  \f'  '*1' 

v^.  V3,  tp) 

df2 

- 

- 

(/»3) 


(5)  Tho  desired  changes  In  6 x (o),  d^,  df!  arc  deter- 
mined from 


“5  X (0)' 
dijj 

■ - e 

Xq  - X (0) 

l 

a 

(^I'O 


since  the  TPBVP  Is  completely  solved  when  the  error 
vector  on  tho  right  side  of  Eq  (kh)  is  driven  to  0, 
or  equivalently,  all  boundary  conditions  are  satis- 
fied, 

(6)  Inverting  tlic  transition  matrix  of  step  4 and  using 
tho  values  of  iS  x (O) , and  dJ2  dotermined  In 
step  5i  Eq  (^3)  is  employed  to  solve  for  (S  x (tp), 
dv,  and  dtp. 

(7)  The  ostimates  In  stop  I arc  nov/  updated  according 
to  the  relationship 


■ ^ ' ' ' 

xp 

xp 

6 x(tp)  x(tp)  dtp 

V 

n 

V 

+ 

dv 

tf 

dtp 

- - 

nev^ 

u top  1 



and  the  procedure  Is  repeated  until  tho  right 
side  of  Eq  (Vl)  Is  .irbitrarlly  small. 


(45) 
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Upon  convergence  of  the  neighboring  extremal  algorithm,  u*(t)  nay 
be  determined  through  Eq  (22).  However,  such  convergence  (s  not  always 
guaranteed.  It  Is  actually  based  on  two  unrelated  factors  which  are 
covered  In  Chapter  IV:  accuracy  of  the  initial  estimate  and  sensitivity 

effects . 

While  both  of  those  factors  v/ere  encountered  to  some  degree,  ax- 
treiiiely  accurate  results  were  obtained  through  a coniblnation  gradient- 
neighboring  extremal  algorithm,  (lore  importantly,  these  results  pro- 
vided the  necessary  basis  for  the  investigation  of  the  non-coplanar  or- 
bital transfer  problem. 
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III-  Non-Coplnnnr  Transfer 


Coord  I na  te  Sys  tem 

In  cl  natural  extension  of  the  coplanar  transfer  problem  derived  in 
the  preceding  chapter,  a spherical  coordinate  system  Is  introduced  as 
shown  in  Fig.  2.  The  coordinate  z is  normal  to  the  Initial  orbital 
plane  with  the  angle  described  in  the  usual  manner.  The  orbits  are 
coapsidal  in  that  a clockwise  rotation  through  the  inclination  angle  I 
about  tile  senil-latus  rectum  of  the  initial  orbit  results  In  a common 
apsidal  linn  for  both  orbits. 

V/hile  tills  restriction  does  not  allow  for  the  most  general  orbital 
transfer  case,  it  Is  imposed  here  as  an  initial  attempt  to  extend  this 
problem  to  out-of-pianc  transfcis.  Moreover,  the  numerico)  technique 
to  bo  presented  could  be  used  to  eventually  handle  the  most  genera i non- 
coplanar  transfer  problem.  , 


Fig.  2.  Non-Coplanar  Coordinate  System 
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Formal ntlon  of  the  Non-Conlannr  TPGVP 


While  this  development  proceeds  In  a similar  manner  to  Chapter  II, 
the  spherical  tr Iqonomctr Ic  relationships  necessary  to  formulate  this 
particular  TPBVP  are  much  more  complex.  For  this  reason,  many  of  the 
details  Involved  In  the  actual  derivation  of  this  problem  are  contained 
in  Appendix  B. 

‘ The  state  equations  are  as  follows 


r « Vr 

(W) 

9 

r 

(^7) 

m 

Vr  - lil  . yii  ,in2  ^ - H- 

(hS) 

r r 

m^  - mt 

Urt  = -VrVO  . 2V4V9  cot  4. 

(50) 

r r 

(m^j  - lilt)  sin  (j» 

-VrViJi  , Vf)^  sin  cos 
V<p  “ + ■ ■ — — ~ — — 

± + ^ 

(51) 

^ niy  - mt 

where  £,j , ^2 1 of'd  are  the  direction  cosines  between  the  thrust  vector 
and  the  r-0-rJ)  coordinates,  respectively.  They  are  Introduced  In  lieu  of 
the  angle  control  of  Cliapter  II  For  ease  In  establishing  optimality  cri- 
teria later, 

In  addition  to  specifying  p^,  e^,  pj:,  e^;,  and  as  In  the  coolanar 
problem,  the  desired  inclination  change  1 must  also  be  selected.  Ini- 
tial conditions  are  determined  by 


r * — 

® I + c_  cos 
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sin  0^ 

Mr  n -<?-  o . o 

o 


(53) 

(5M 


ve.  - 


(55) 


- 0 


(56) 


with  terminal  constraints 


where 


Pf 


1 + Cf  cos 


- 0 


f 


ii/g  ■ (j)f  - [()q  + arctan  (tan  I cos  0jr)l  * 0 
f 


Cehr  sin  Q'f 

a MVf  - -i— ' f.  » 0 


Pf 

» V0^  - ^ (cos  I - sin  1 cos  9^  cot  dj.)  » 0 
hf  sin  I sin  6f 

|J>5  « V^|T  + 0 


sin  (arctan  (tan  I cos  Of)) 

O'c  » arccos  , I 

^ . sin  i i 


Again  the  Mayor  formulation  of  the  cost  function  gives 


J “ ti 


(57) 

(58) 

(59) 

(60) 
(61) 

(62) 


(63) 


and  the  Hamiltonian  becomes 
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H » X Vr  4-  X -—  + ,\ 

I 2 r 3 r 


+ X ( V:^  + yq.?_s_I nij.  _ ^ ]h_  \ 

^ f ••  r2  m,.  - mt  ' 

f-VrVe  2VjV0 


+ X 


cot  6 


U. 


(a,.  - mt)  sin  (Ji 


] 


+ X, 


-i- ) 

-ml*  ' 


r.iQ  - mt 


This  leads  to  the  costate  equations 

. v|  X * n.  >.,*  . 2j.N  X. 

* 2 r2  ' \r2  rV  ‘ 


/•VrVJ  2V3V:  cot  cS\  , 

(■~T'  --7“-^)  ^5 


/“V_r\,'j:_  + V '1^'  sin  ^ COS  ;i>  \ 

V'  ^ — •— Y' 


''c 


x„  - 0 


~2  j l_n  _'._co_L:LJ/JL1  1 
r 4 


r;_2y?yj)_ 

l_r  sin2 


Tij  cos  ij> 


4*  — — ^ 

("1^  - nil)  sin^  -!>J 


■|.  _(^rJn^  iji  “ cos^  <{i)  ^ 


; . V + '/O  . . Vi  > 

^4  'h  ^ r “ '^5  r” 


" " >'  '^2 


1 . 2Vt)  sin^  .i 


\ * Cr* 


•J.\l  0_  s j_n_  y_ _4.  ^ 


ftoproduced  ,from 
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6 r3r4  r srg 


(70) 


and  transversal i ty  cond I tions 

Xj(tf)  “ + v^  (cos  I “ sin  1 cos  Of  cot  i{if) 


Vs  sin  I sin  0|! 


(71) 


^2  (^ f ^ 


Pf^f  sin  0 ^f  38'*^ 

^ (1  + e^;  cos  0'|;)^  30^ 


sin  e,  tan  I 

+ V,  £ 


- V. 


cos  Q'f  38'p 


I + tan^  1 cos^  9^ 

- Vi,  — f-  sin  1 sin  0,  cot  (J>. 

- tf  f f 

hf 

+'  Ve  --  sin  t cos  Q( 

hf  sin  1 cos  0f 

A,(tf)  ■ V2  “ \>4 pT 

3 ^ rf  sln2 


38. 


(72) 

(73) 


X4  (tf ) ■ V3 


■\5(tf)  ^ V4 


Xg(tf)  - Vg 


H(tjr)  ’ “ 1 


(74) 

(75) 

(76) 

(77) 


Minimizing  the  Hamiltonian  with  respect  to  the  controls  li,  ^2* 
and  i^3  yields 


i.  * 


/X^2  4.  (Aj/sin  <>)^  + X(j^- 


(7B) 
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® XI 


stn  Xj  cos  Xj 


(mIq  - mt)  A 


Xg  Xg 

^ ‘ ~T  ^2  ^ 7T  S 

Xj  x^ 


( 


(■ 


Xg^  + Xij^  sin^  X3  2p 

Xj2 

-X4X5  - 2X3X5  cot  X3 


r) 


ij  ^ 


) 


-x^Xg  + X;j2  sin  X3  cos  Xj,^  ^ 

+ [ — ' ->■•  '■  ■■  i.i 1. 1-  III  "-  ■ j \ 


< 


'I 


Xj  - 0 


sin  X-  cos  x.x 


,1'^S 


•2xcX 


5^6 


Xj  sin^  X3 
TXg  cos  X3 

(nip  - iTit)  A sln^  Xj  J 


^5  ' 


X6(0)  - 0 


(87) 


hf 

X^(tf)  - Vj  f (cos  1 

^If 


-sin  I cos  K^f  cot  x,j;) 


sin  I sin  x^f 
^If 


+ 


.)  - 

■ 

s i n 

1 + 

tan' 

hf 

s 1 n 

^Tf* 

5 1 n 

(83) 

5 1 n 0 SO 


I cos  X 


2f 


(89) 


X3(tf)  - Vj 


- Vk 


lip  sin  I cos  Xjp 

Xjf  Sln^  Xj^r 


<tproaucBa  from 
3c»t  «v«ll»ol*  toRV: 
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Xj^CsIn^  Xj  - cos^  Xj). 
^1 

^4  - - h ^ r ^5  + "■ 

Xj  X^ 


1 


X - 


2x5  sln^  xj 


5 X2  2 

/X4  2X(3  cot  X3  V 

— ) , 

2x5  sin  X3  cos  X3 
> X, 

• 1 2Xg 

^6  " “ Xj"  “ Xj" 


. X 

Xi  4 


2xg  cot  X3 


‘1 


><4 

^5  ^5 


X,,(tf)  - V3 


■ \ 


X^ltf)  - V5 


♦1  ■ *if  ■ i~~f“c5r»7  ■ “ 

'^2  “ ^3f  ” ^'*’0  ^ (tan  I cos  x^p)]  ■ 0 

'K  " 


(90) 

(91) 


(92) 


(9:) 

(9A) 

(95) 


epli|  sin 
%f  " *P|T“ 

!v 

■ 

0 

(96) 

X c - (cos 

5f  xif 

1 “ sin  1 cos  Xgp  cot  Xjp)  “ 0 

(97) 

X . ' sin 

6P  Xip 

1 sin  x^p  « 0 

(9«) 

n - H(tp)  + I “ 0 


(os) 
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12  constants  of  Integration 
+5  unknown  multipliers  v, 

til 

« lO  unknov^ns 


6 Initial  conditions 
+5  terninal  constraints 
+ / transvnrsa 1 1 ty  conditions 

= 13  equations 


As  in  the  coplanar  transfer  TPiiVP,  sufficient  relationships  exist 
to  obtain  numerical  solutions  for  the  non-coplanar  problem.  Further- 
more, it  is  apparent  that  a snecified  Inclination  change  of  0 reduces 
this  TPBVP  to  the  coplanar  case. 

No  I nhbor Ing  Extremal  Algorithm 

This  similarity  suggests  the  likelihood  of  obtaining  solutions  to 
the  non-coplanar  TPBVP  by  Initially  raking  a small  Inclination  change 
to  an  optimal  coplanar  transfer  trajectory.  The  use  of  a grocilont  al- 
gorithm to  obtain  a reference  Input  for  the  neighboring  extremal  al- 
gorithm can  thus  be  obviated. 

The  particular  neighboring  extremal  algorithm  used  is  identical 
In  nature  to  that  outlined  In  Chapter  II  and  for  that  reason  will  not 
be  detailed  again  here.  In  fact,  the  only  changes  arise  through  an 
Inclusion  of  the  four  added  initial  and  tcrninal  conditions  to  this 
enlarged  12  state  system. 

The  original  technique  of  specifying  a small  Inclination  change 
with  on  optimal  coplanar  trajectory  as  a reference  Input,  proved  quite 
successful  in  obtaining  accurate  numerical  solutions  through  this 
neighboring  extremal  algorithm.  Moreover,  once  a solution  is  deter- 
mined for  a small  inclination  change,  it  Is  then  possible  to  use  this 

solution  as  a reference  Input  and  Increase  the  soutlfled  Inclination 

» 

change.  This  procedure  Is  repeated  to  obtain  numerlcol  results  for 
relatively  large  values  of  inclination  change. 
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The  major  drawback  of  this  approach,  however,  Is  the  excessive 
■number  of  iterations  required  for  convcrijence,  especially  In  view  of 
sensitivity  effects  which  may  severely  limit  the  size  of  the  variation 
in  inclination.  This  problem  will  be  discussed  in  the  next  chapter, 
where  a curvo-f  I tt  ln;j  technique  is  presented  and  shown  to  be  very  ef- 
fective in  ovcrconing  this  obstacle. 
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C ompu  t e r I nip ) cnicn  t_a  1 1 on 

In  applying  the  numerical  algorithm  of  Chapters  II  and  III,  two 
computer  programs  were  written:  CPTRAN  for  the  coplonar  transfer  and 

NCPTRAH  for  the  non-coplanar  problem.  A separate  listing  for  those 
two  programs  Is  contained  In  Appendix  C and  D,  respectively. 

For  practical  reasons  alone,  the  seml-latus  rectum  of  both  the 
Initial  and  final  orbits  remains  fixed  In  this  analysis.  V/lth  one 
exception,  transfer  from  the  Initial  orbit  was  always  specified  to 
stort  at  perigee.  Thus  the  only  parameters  that  varied  were  the  eccen** 
tricitics  of  the  Initial  and  terminal  orbits  as  well  as  the  Inclina- 
tion angle  In  the  non-coplanar  problera. 

Geocentric  canonical  units  were  used  for  all  computational  work 
with  the  results  displayed  accordingly.  Thus,  the  unit  of  length  was 
the  radius  of  the  earth  and  the  gravitational  constant  u was  normalized 
to  one.  This  resulted  in  a time  unit  (TU)  of  13*^i*<7  minutes.  Vehicle 
mass  was  also  normalized  to  one  to  produce  a constant  mass  flow  rate  of 
.005.  In  all  cases,  accuracy  was  required  to  the  fifth  decimal  place 
before  convergence  of  a solution  was  presumed. 

An«alysis  of  the  Copl.an.'ir  Problem 

As  previously  mentioned,  the  effectiveness  of  the  no Ighbor Ing  ex- 
tremal algorithm  Is  based,  to  a large  extent,  on  the  reference  Input 
provided  through  the  gradient  method.  It  was  thus  desirable  to  achieve 
approximate  solutions  to  the  coplanar  TPDVP  from  the  gradient  algorithm 
in  as  few  Iterations  as  possible.  This  was  accomplished  In  the ‘follow- 
ing manner. 
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The  first  few  Iterations  show  the  auqnicntcd  cost  function  of  Eq 
(35)  to  decrease  rapidly,  after  which  a plateau  Is  reached  where 
further  reductions  In  the  cost  beconc  quite  small.  It  Is  at  this  point 
that  the  b multipliers  of  Eq  (36)  arc  exploited  to  overcome  this  dif- 
ficulty and  to  speed  up  the  convergence  jioccss.  In  CPTRAN  a c^'eck 
was  thus  made  on  the  not  decrease  In  If  this  decrease  was  less  than 

10"'  the  time  consiiinlng  process  of  updating  the  control,  and 

starting  a new  Iteration  was  bypassed.  Instead,  the  b|  were  undated 
accordlnq  to  step  6 of  the  alqorlthi  until  a reasonoble  decrease  In 
Jg  was  achieved.  This  small  adjustment  was  Instrumental  In  obtaining 
results  accurate  to  the  second  decimal  place  In  50  Iterations. 

In  transitioning  to  the  neighboring  extremal  algorithm,  the  final 
values  of  u and  tf  obtained  above  are  used  to  compute  n reference 
trajectory  for  the  coplanar  TPBVP.  The  eight  parameters  needed  in 
step  1 of  the  procedure  are  thus  determined,  In  ovary  case  this  ini- 
tial estimate  was  sufficiently  close  to  the  optimal  for  convergence 
of  the  algorithm. 

The  linear  assumption  inherent  In  Eq  (143)  places  groat  emphasis 
on  the  selection  of  e In  Eq  (l|4).  As  the  boundary  conditions  become 
closer  to  being  satisfied  this  I'inoority  assumption  Is  more  valid  and 
the  value  of  c con  be  Increased  accordingly,  up  to  a ino.ximum  value  of  1. 
On  the  other  hand,  a poor  initial  estimate  completely  destroys  such  a 
supposition,  and  no  matter  how  small  c Is  picked,  uss  of  this  method 
v/ould  be  questionable.  The  value  of  c used  In  CPTRAII  was  thus  con- 
tinuously upd.itcd  by  dividing  .005  by  tlie  root-meum-squaro  of  the 
error  vector  In  Eq  (Vt)  until  e reached  Its  maximum  value  of  1. 
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V/hllc  n complete  sensitivity  analysis  Is  beyond  the  scope  of  this 
‘thesis,  a short  note  on  sensitivity  effects  Is  In  order.  Due  to  the 
nature  of  the  costate  or  Eul er-Lagrange  equations,  as  x increases 
along  on  extremal  path  \ decreases  in  magnitude.  This  property  in  turn 

i 

causes  largo  inacturacles  In  the  trans I tion  matrix  of  the  neighboring 
extremal  algorithm  througli  round-off  and  integration  errors  (Ref  I: 
215).  It  is  then  possible  for  an  overcorrcctlon  In  the  terminal  con- 
ditions of  step  (j  In  the  neighboring  extremal  algorithm  and  divergence 
from  the  solution. 

It  was  found  that  a perturbation  on  tlio  order  of  10"^  In  stop  k 
of  the  neighboring  extremal  algorithm  produced  0 transition  matrix 
equivalent  to  one  produced  using  any  smaller  perturbation,  For  this 
reason  a perturbation  of  .01  was  used  In  CPTRAii  to  generate  tlio  transi- 
tion matrix.  The  sllglit  delay  encountered  In  final  convergence  of  the 
solution  Is  thus  explained  by  the  effect  of  the  above  noted  sensitivity 
In  forming  a transition  matrix  of  sufficient  accuracy  to  improve  the 
previous  estimate. 

The  effectiveness  of  the  combination  grndicnt-ncigliboring  extre- 
mal algoritlim  is  that  no  knowledge  of  the  solution  Is  regulrod  before- 
hand. For  any  Initial  estimate  of  the  control,  no  matter  how  poor, 
accurate  solutions  to  the  TPDVP  are  ultimately  achieved.  However, 
this  flexibility  resulted  in  a penalty  of  Increased  computer  time. 

Thus,  once  an  optimal  solution  was  obtained  for  a particular  sot  of 
orbital  parameters,  it  was  more  efficient  to  use  this  solution  as  an 
estimate  In  the  neighboring  extremal  algorl thm  wfion  small  variations 
in  these  orbital  parameters  were  desired.  In  some  Instances,  such  as 
a small  cliOnge  In  the  eccentricity  of  the  terminal  orbit,  It  was  then 
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possible  to  bypass  the  gradient  portion  of  this  method.  This  caused 
a 50i5  reduction  in  the  computational  time  required  for  convergence. 

Ites u 1 ts  for  _t j]e.  Coplanar  Pro bljBm 

Transfer  between  two  circular  orbits  v;as  chosen  as  an  obvious 
starting  place.  The  well  known  spiral  transfer  trajectory  was  obtained 
as  shown  In  Fig.  3i  with  the  control  history  plotted  as  a function  of 
transfer  time  In  Fig.  A.  This  steering  program  compared  favorably  vylth 
that  obtained  by  Faulders  in  a similar  problem  (Ref  2:A5).  On  all 
figures  PO,  EO,  PF,  and  EF  refer  to  the  semi-latus  rectum  and  occen- 
trlcity  of  the  Initial  and  terminal  orbit,  respectively. 

Results  for  various  combinations  of  eccentricity  were  also  cb- 
talnod,  As  one  might  suspect,  there  is  a remarkable  similarity  in  con- 
trol histories  when  eccentricities  are  the  same  or  even  approximately 
so.  This  effect  Is  demonstrated  by  comparing  Fig.  A witli  the  control 
lilstory  plots  of  Figs.  5 and  (i.  Wote  that  all  three  cases  arc  charac- 
terized by  a rapid  reversal  in  thrust  direction  through  the  halfway 
point  in  transfer  time. 

This  is  In  sliarp  contrast  to  orbits  of  differing  eccentricity. 
Corresponding  trajectories  and  steering  programs  Tor  two  such  cases 
are  displayed  In  Figs.  7 thru  10.  Fig.  8 shows  a delay  in  the  thrust 
reversal  noted  above  whereas  tlio  opposite  cFFoct  Is  observed  In  Fig. 

10.  That  Is  to  say,  In  transfering  to  an  orbit  of  greater  eccentri- 
city the  net  effect  Is  a shift  of  the  control  plot  to  the  right. 
Similarly,  a shift  to  the  left  results  wlien  the  Initial  orbit  Is  of 
larger  eccentricity. 

Fl«,.  9 shows  the  only  case  considered  whore  departure  occurred 
from  other  than  perigee,  Oq  “ A5“  to  bo  exact.  In  this  particular 
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Instance  a minimal  change  In  the  control  law  resulted.  However,  total 
time  of  transfer  was  some  .71  TU's  less  than  the  same  problem  v/here 
departure  was  specified  at  perigee.  This  fact  Is  mentioned  here  only 
to  show  that  no  general  statement  can  be  made  concerning  transfer  time. 
It  Is  strictly  0 function  of  the  eccentricities  of  the  orbits  and  the 
starting  point  from  the  Initial  orbit. 

Analysis  of  the  Hon“Coplonar  Problem 

In  Chapter  III  the  recursive  method  of  obtaining  solutions  for 
Increasing  values  of  Inclination  between  specified  orbits  was  exolelned. 
This  technique  v/as  Initially  employed  here  Inasmuch  as  reference  solu- 
tions were  already  available  fron  the  coplanar  analysis.  Very  slow 
convergence  of  a gradient  algorithm  In  this  twelve  state  system  with 
three  control  variables  was  olso  anticipated. 

On  the  other  hand,  solutions  for  an  Inclination  change  approach- 
ing 90"  wore  an  ultimate  objective.  An  alternative  approach  was  thus 
necessary  to  avoid  the  arbitrariness  of  the  recursive  technique.  For 
Instance , after  successfully  achieving  an  Inclination  change  (IHC)  of 
36“,  this  procedure  would  only  work  for  an  Increase  In  Inclination  of 
loss  than  ono-fourth  of  a degree  using  the  3^*  transfer  as  the  Initial 
estimate.  This  was  clearly  unacceptable. 

The  approach  that  was  adopted  v*(as  stralqhtforv/ard  and  very  ef- 
fective In  surnrauntlng  such  uncertainty.  It  consisted  of  using  the 
twelve  parameters  Xf,  v,  ond  tf  obtained  from  tlie  previous  three  solu- 
tions for  respective  values  of  Inclination,  A least  square  curve-fIt 
was  then  made  to  estimate  these  par.smcters  for  the  now  value  of  In- 
clination. To  Illustrate,  consider  the  above  nentlonod  problem  where 
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the  recursive  technique  was  stalled  at  36**  The  solutions  obtained 
for  inclination  change  of  30* » 33*i  and  3<j*  provides  reference  points 
for  each  of  the  twelve  parameters  needed  In  stop  1 of  the  nolrihboring 
extremal  algorithm.  Curve-fitting  thus  provided  a second  order  poly- 
nomial which  In  turn  produced  an  estimate  of  x*p,  v",  and  t^  for  an  in- 
clination change  of  39"* 

Figs.  11  and  12  show  a plot  of  these  parameters  as  a function  of 
Inclination  change  for  a transfer  between  the  same  circular  orbits 
used  In  the  coplanar  analysis,  uesldes  displaying  the  obvious  non- 
linear cl^aracter I St  Ics  of  the  problem,  these  graphs  may  be  Interpreted 
as  sensitivity  plots,  since  they  show  the  effect  of  a sr.al)  variation 
in  Inclination.  Note  that  while  small  changes  In  Inclination  resulted 
in  relatively  snail  changes  In  and  t^,  such  Is  not  the  case  for  *v. 
Again  using  the  3(j’  value  of  Inclination,  It  can  bo  observed  from  Fig. 
12  that  V21  Vj,  Vt,,  and  Vg  are  rapidly  changing  for  even  a small  In- 
crease In  Inclination.  This  explains  the  difficulty  encountered  In  the 
recursive  method  since  results  for  an  Inclination  change  of  36*  are 
thus  shov/n  to  provide  a poor  estimate  for  .my  lornor  value  of  Inclina- 
tion. 

As  expected,  upon  approaching  a specified  Inclination  change  of 
35"  convergence  became  extremely  slow.  This  was  so  for  two  reasons: 
the  extremely  long  transfer  times  caused  roun<l-off  and  Integration 
errors  to  have  a much  greater  effect  than  In  the  coplanar  problem 
where  this  factor  was  discussed,  and,  at  an  Inclination  change  of  90* 1 
the  state  equations  go  umloftned  as  is  appronchos  0.  This  In  Itself 
causes  extreme  sensitivity  to  any  varl.atlon  In  0 near  optimal  trajec- 
tory, The  overall  result  was  a complete  breakdown  In  the  neighboring 
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Variation  of  Xf  to  Inclination  Change  for  Transfer  Between  Circular  Orbits 


Fig.  12.  Variation  of  v and  t,  to  Inclination  Change  for  Transfer  Between  Circular  Orbits 
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extremal  algorithm  and  for  this  reason  no  solutions  were  obtained  be- 
yond an  inclination  change  of  S5*.  Actually,,  this  apparent  limitation 
could  have  been  overcome  by  a simple  rotation  of  the  coordinate  system 
to  avoid  this  singularity. 

I 

Results  for  the  Non-Copl anar  Probl em 

For  ease  In  understanding,  the  control  histories  are  displayed 
graphically  as  tv.o  tlirust  directions.  Ul,  In  keeping  with  the  co- 
planar  problem  for  comoarative  purposes,  represents  the  thrust  direc- 
tion as  measured  clockwise  from  tlic  local  horizontal  In  the  r-0  plane. 
U2  Is  then  defined  as  the  thrust  direction  measured  clockwise  out  of 
the  r-0  plane.  Previous  notation  Is  used  on  all  figures  with  the 
addition  of  IIIC  to  denote  Inclination  change  In  degrees. 

Fig.  13  shows  the  control  histories  for  a snail  out-of-plano 
transfer  between  tl)c  same  circular  orbits  used  In  the  coplanar  analy- 
sis. llotc  that  Ul  Is  almost  Identical  to  tlie  control  shown  In  Fig.  h, 
as  should  be  the  case.  U2  Is  also  reasonable  since  very  little  out-of- 
plane  thrusting  would  be  required. 

As  demonstrated  in  Figs.  and  15,  a significant  change  In  con- 
trol histories  occurred  as  Inclination  was  increased.  In  each  case, 

Ul  bears  little  rosonblance  to  any  steering  program  obtained  In  the 
coplanar  problem,  since  most  of  the  thrust  Is  now  directed  out-of- 
plane.  On  the  other  hand,  U2  displayed  a remarkable  similarity  for 
an  inclinotlon  change  of  30*,  even  thougli  tlie  transfer  was  between 
circular  orbits  In  Fig,  1 ^ and  widely  varying  elliptical  orb  I ts  In 
Fig.  15. 

Those  results  suggest  the  follovilng  relationships:  out-of-plane 

control  IWstory  (U2)  is  primarily  a function  of  specified  Inclination 
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change  while  the  In-plane  steering  program  (U1)  varies  v/Ith  a change 
In  eccentricity  in  much  the  same  manner  as  a coplanar  transfer.  Fig, 

16  Is  further  confirmation  of  these  relationships. 

Further  Inspection  of  Figs.  lA-16  reveals  tliot  U2  is  not  unlike 
a bang-bang  controller  in  that  the  out-of-plane  thrust  Is,  for  the 
most  part,  nearly  normal  (±  00°)  to  the  instantaneous  orbital  plane  of 
the  transfer  trajectory.  This  particular  result  was  obtained  In  all 
transfers  of  3O°-A0“  Inclination  change  and  Is  not  surortslng  In  view 
of  the  coordinate  system  used,  it  is  Mentioned  here  because  transfers 
In  tiiis  Inclination  range  are  often  required  to  achieve  a polar  or 
equatorial  orbit  for  vehicles  launched  from  the  United  States. 

Graphical  presentation  of  transfer  trajectories  was  prohibitive 
in  view  of  the  added  dimensionality  of  the  non-coplanar  problem.  How- 
ever, for  values  of  Inclination  greater  than  transfer  trajectories 

in  excess  of  one  earth  revolution  were  obtained.  Fig,  17  portrays  a 
typical  control  history  for  a transfer  of  this  type.  Due  to  the  long 
transfer  time  required,  the  previous  bang-bang  property  of  U2  Is  no 
longer  present. 

Unlike  the  coplanar  problem,  transfer  time  between  non-copl anar 
orbits  is  very  much  dependent  on  required  inclination  change.  In  fact, 
the  TF  plot  of  Fig.  12  shows  transfer  time  as  a.  monoton  lea) ly  Increas- 
ing function  of  specified  Inclination  change  for  the  given  circular 
orbi ts. 
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The  problem  of  minimum  time  transfer  to  some  desired  terminal 
orbit  using  a constant  magnitude,  low  thrust  propulsion  system  has 
been  Investigated. 

For  the  case  of  a conlanar  transfer,  the  application  of  a combina- 
tion gradlcnt-nclgtiborlnq  extremal  algorithm  was  shown  to  be  a most 
effective  moans  of  obtaining  numerical  solutions  to  the  resulting  non- 
linear TPBVP.  The  particular  appeal  of  tills  method  was  that  no  know- 
ledge of  the  solution  was  required  a.  to  obtain  extremely  ac- 

curate results  Thus,  a combination  gradient-neighboring  extremal 
nhjortthn  would  seem  to  be  applicable  to  a wide  class  of  orbital 
transfer  problems. 

On  the  other  hand,  for  greater  efficiency  In  computer  time,  It  Is 
possible  to  bypass  the  gradient  portion  of  tills  method  once  an  optimal 
solution  Is  obtained  for  transfer  betv/oen  specified  orbits.  This  Is 
accomplished  by  using  the  above  solution  os  o direct  input  to  the 
neighboring  extremal  algorithm  when  making  some  small  variation,  such 
as  Increasing  eccentricity.  In  the  given  orbital  parometers. 

V/lth  this  opprooch,  the  coplanar  onalysts  was  Instrumental  In 
extending  the  problem  to  non-coplanar  terminal  orbits.  Upon  achiev- 
ing solutions  to  the*,  non-coplanar  problem  for  small  values  of  Inclina- 
tion, a method  of  using  those  Initial  results  v/as  developed.  It  con- 
sisted of  generating  now  cstihiatcs  to  the  neighboring  extremal  algo- 
rithm by  a polynoml<il  approximation  of  poramoters  as  n function  of 
Inclination.  In  this  manner,  solutions  to  the  non-coplanar  problem 
were  achieved  bo  tween  orbits  of  widely  varying  I ncl  ln.it  Ion. 

h7 
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In  view  of  the  accuracy  of  these  results,  this  approach  forms  the 
basis  of  an  effective  numerical  technique  for  completely  general  opti- 
mal low  thrust  orbital  transfer. 
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Derivation  of  Equations  for  Coplanar  Tronsfer 
State  liquations 

Usln<]  the  standard  vector  notation  of  o’^  and  e,j  to  denote  unit 
vectors  in  the  r“0  plane,  the  position  vector  of  the  vehicle  at  any 
time  may  be  written  as 


~ - rT^  (\“l ) 

wlLl.  the  anquiar  velocity  of  the  coordinate  system  given  by 

u - 9e'n  (A-2) 

where  e„  is  a unit  vector  normal  to  the  r-0  plane.  Differentiating 
£q  (A-I)  twice  yields 

« 

r - ro~r  + rOOg  (A-3) 

r m (r  - rO^)  ‘c^  + {z'ri  + ro)  tTg  (A-h) 

as  vehicle  velocity  and  acceleration,  respectively. 

The  restricted  two-body  cquacio.-i  v.'lth  terms  to  account  for 
vehicle  thrust  has  the  form 


«« 

'r 


/ ly.  + ! sin  u'^  Op  + cos  u Cg 

V hIq  “ mt  y niQ  - mt 

(A-5) 


where  u is  measured  clockwise  from  the  Cg  direction. 

Now  lot 

r ■ Vr  (A-6) 
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and  equatinq  vector  components  on  the  right  hand  side  of  Eqs  (A-M  and 
(A-5) 


Vr  ^ sin  u 

**  r“  - mt 

(A-8) 

• -VrVO  ^ T 

VO  ■ + cos  u 

' '"o  “ 

(A-9) 

which  arc  the  four  state  equations. 

Doundary  Conditions 

The  general  conic  solution  to  the  restricted 

two-body  equation 

has  the  form 

r R 

1 + e cos  0 

(A-10) 

Differentiating  Eq  (A-IO)  y’olds 

Vr  « rJi-Siili. 

P 

(A-I 1) 

whore  h is  conserved  specific  angular  momentum  such  that 

h • r^  0 

(A-12) 

and  applying  Eq  (A-7)  In  Eq  (A-12)  results  In 

VO  - - 

f* 

(A-13) 

Thus,  Eqs  (A-IO),  (A-II),  and  (A-13)  are  tho  three  conditions 


which  must  be  satisfied  In  both  the  Initial  and  terminal  orbits. 
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Derivation  of  Egua^t Ions,  for  Non-Coplanar  Transf ej 


State  Equations 

The  standard  vector  notation  of  o'^,  oq,  and  Is  atialn  used  to 
denote  unit  vectors  In  the  r-0“<{i  spherical  coordinate  system,  which  has 
angular  velocity 


u)  ■ 0 cos  cq  - 6 sin  1^ 

D 1 fforont latlnq  the  position  vector,  r^j.,  gives 

A ^ ^ 

r ■ re^  + rO  sin  Gq  + rij  e,^ 

us  vehicle  velocity  and 
« 

• t ^ 

r ■ (r  - r«j2  - rO^  sln^4)  o’,. 

(rO  sin  (Ji  + 2r0  sin  <j»  + Zrifo  cos  ij»)  o’q 
+ (r^  + 2rij  - rO^  sin  ^ cos  <}•)  (H"3) 


(B-1) 


(B-2) 


as  vehicle  acceleration. 

The  restricted  two-body  equation  with  additional  terms  for 
vehicle  thrust  becomes 


r 


Vr2 


Tt, 


m^j  - mt 


) 


n,  _ 

+ i--  c„ 

- mt 


+ 


iHq  “ mt 


(o-»0 


t 

where  tp  tj,  and  <.3  arc  the  direction  cosines  between  the  thrust 
vector  and  e,.,  o'^ , and  e^,  respectively. 
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In  addition  to 

r « Vr  (B-5) 


lot 


(B-6) 


(B-7) 


and  again  oquating  acceleration  terms  in  Eqs  (B-3)  and  (D“^0 
Vr  - s|„2^  - iU  ^ 


nirt  " mt 


VO . :ysiil . + li.?. 

^ •'  (it)_  “ lilt)  sin  ij» 


v^  - ±ry±  ^ 'i9.LjJ]L±sp±J.  + 

•*  f"  niQ  - mt 


which  aro  tho  six  state  equations. 


(B-8) 

(D-9) 

(B-10) 


jiLli J ai  Conditions 

Three  of  tlic  five  unspecified  initial  conditions  are  the  same  as 
those  derived  In  Appendix  A,  Eqs  (A-IO),  (A-II),  and  (A-13)*  tp  of  tho 
Initial  orbit  Is  fixed  at  90®  In  view  of  tho  definition  of  tlio  coor- 
dinate system  In  Chapter  III.  This  In  turn  sots  tlio  value  of  Vij)  at  0 
cverywhwro  In  tho  initial  orbit. 


Gond  1 1 Ions 

The  establishment  of  terminal  conditions  for  an  Inclined  ellipti- 
cal orbit  requires  tho  lisk  of  spliorical  tr Iqonomotric  rcl at  lonslilps 
to  solve  for  tho  angles  i(i'‘  and  0'  ns  shov/n  In  Fiq.  10. 

From  ilapior's  ruio  for  right  splioric.il  trinnolos 
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The  final  two  conditions  are  established  through  the  definition 
of  conserved  specific  angular  momentum 

h *•  r X r “ rV<^  cq  - rVO  sin  (}ie^  (D-18) 

In  the  terminal  orbit i h may  be  written  as 

h « h sin  I ej(  + h cos  1 Oj,  (B“19) 

V'/hcro 

e*^  “ cos  0 sin  (|i  Op  + cos  0 cos  4>  - sin  0 Cq  (B-20) 

■ cos  (j)  ep  - sin  ^ (B“21) 

After  substitution  of  Eqs  (El-20)  and  (B-21),  Eq  (0-15)  Is  equated 
to  Gq  (Q-lQ)  to  yield 

VO  ■ ~ ( cos  I - sin  1 cos  0 cot  i^)  (0-22) 

V(|)  • p-^‘  sin  I sin  0 (B-23) 
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Appendix  C 

The  following  contains  a listing  of  the  computer  program,  CPTRAN, 
which  uas  used  to  obtain  results  to  the  coplanar  TPBVP  of  Chapter  II. 
Bes ides • those  subroutines  listed,  SET,  STEP,  and  MINV  were  also  em- 
ployed. SET  and  STEP  were  called  for  integration  purposes,  wfille  MINV 
was  needed  to  Invert  the  transition  matrix  of  the  nelgliboring  extremal 
algorithm.  All  three  of  these  subprograms  are  found  on  the  general 
package  of  AFITSUBROUTINES  which  Is  on  permanent  file  at  the  Air  Force 
Institute  of  feclinology. 
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Appendix  0 


The  following  contains  a listing  of  the  computer  program,  tICPTRAN, 


which  was  used  to  obtain  results  to  the  non-coplanar  TPBVP  of  Chapter 
III.  The  AFITSUCROUTINES,  SET,  STEP,  and  MI.'IV  were  again  employed  In 
the  same  manner  as  explained  In  Appendix  C. 
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Vita 

Eugene  A.  Smith  v(as  born  27  July  UV»5,  in  Micldlovlllo  NY.  He 
attended  West  Canada  Volley  Control  School  where  he  was  valedictorian 
of  his  19*j3  graduating  class.  An  Air  Force  Academy  appointment  fol  ■ 

t 

lowed,  and  upon  graduation  In  lOIj?.  ho  was  commissioned  a second 
lieutenant  and  sent  to  Reese  AFD  for  Undergraduate  FMlot  Training. 
After  UPT,  he  served  In  South  East  Asia  flying  A-l's  and  0-2's.  From 
this  overseas  tour  he  received  the  Silver  Star,  the  Distinguished 
Flying  Cross,  and  cinht  Air  Medals. 

Upon  returning  to  the  United  States,  Lt  Smith  was  assigned  to 
VMllIams  AFC  as  a T-.32  Instructor  pilot.  H«  Instructed  for  two  years 
and  was  In  Cfieck  Section  for  one  year. 

In  August  of  l'-l72  Captain  Smith  was  nottFled  of  acceptance  to 
tl>e  Air  Force  Institute  of  Technology  In  the  Astronaut  leal  Engineering 
progr.am. 

Captain  Si.iltli  curreiU  ly  lives  In  Forest  Ridge,  Dayton  OH  with  his 
wife,  Glenda  and  daunlitcr,  Andrea. 


This  thesis  v/ns  typed  hy  Katherine  P.in'Jall 
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